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ABSTRACT 

We take up from a library of 12 galaxy clusters featuring extended X-ray observations of their Intra 
Cluster Plasma (ICP), analyzed with our entropy-based Supermodel. Its few intrinsic parameters — 
basically, the central level and the outer slope of the entropy profile — enable us to uniformly derive 
not only robust snapshots of the ICP thermal state, but also the 'concentration' parameter marking 
the age of the host dark matter halo. We test these profiles for consistency with numerical simulations 
and observations. We find the central and the outer entropy to correlate, so that these clusters split 
into two main classes defined on the basis of low (LE) or high (HE) entropy conditions prevailing 
throughout the ICP. We also find inverse correlations between the central/outer entropy and the halo 
concentration. We interpret these in terms of mapping the ICP progress on timescales around 5 Gyr 
toward higher concentrations, under the drive of the dark matter halo development. The progress 
proceeds from HEs to LEs, toward states of deeper entropy erosion by radiative cooling in the inner 
regions, and of decreasing outer entropy production as the accretion peters out. We propose these 
radial and time features to constitute a cluster Grand Design, that we use here to derive a number 
of predictions. For HE clusters we predict sustained outer temperature profiles. For LEs we expect 
the outer entropy ramp to bend over, hence the temperature decline to steepen at low z; this feature 
goes together with an increasing turbulent support, a condition that can be directly probed with 

the SZ effect. We finally discuss the looming out of two intermediate subsets: wiggled HE at low z 
that feature central temperature profiles retaining imprints of entropy discharged by AGNs or deep 
mergers; high-z LEs, where the cosmogony /cosmology had little time to enforce a sharp outer entropy 
bending. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (A399, A1656, A2218, 
A2256) — X rays: galaxies: clusters — method: analytic. 
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1. INTRODUCTION 

Rich galaxy clusters, with their overall masses M ~ 
10^"'^ Mq and large virial sizes R ~ Mpcs^, constitute the 
most recent cosmic structures with high contrast, but 
still developing at low redshifts. Their gravitationally 
dominant dark matter (DM) halos contain an apprecia- 
ble amount m « 0.16 M of hot, diffuse baryons in the 
form of an intracluster plasma (ICP) at virial tempera- 
tures ~ G Mrup/ 10 R ~ several keVs and with av- 
erage densities n ~ 10"'^ cm~'^. The ICP conditions can 
be probed in X-rays through its strong bremsstrahlung 
emissions of powers Lx oc r? T^l"^ R"^ « 2^q44-45 gj.g g-i^ 

Our main focus here will be on the physics of the ICP, 
and specifically on its 'entropy' 



k = ksT/n' 
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or better adiabat (see Bower 1997), which is simply re- 
lated to the true specific entropy s by As = 3/2 In fc. The 
quantity k will conveniently constitute our leading state 
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variable, due to its basic properties: it is eroded and even- 
tually erased at the cluster center by radiative cooling; 
it is produced at shock fronts driven both by supersonic 
infiows across the cluster boundary and by central out- 
flows; it is conserved and stratified upon adiabatic com- 
pression of the outer intergalactic medium (IGM) into 
the ICP contained by the DM potential well. 

Our scope will be to relate the entropy levels in the 
ICP to the evolution of the containing DM halos. As 
to the latter, we will refer to the standard scenario in- 
cluding hierarchical formation and secondary infall, up- 
dated by state-of-the-art iV— body simulations and an- 
alytical works (e.g., Zhao et al. 2003, Fakhouri et al. 
2010, Wang et al. 2011). This scenario envisages a first 
stage of fast collapse and major mergers forming the halo 
bulk from the top of the initial density perturbation; this 
is followed by a slow development of the outskirts by 
accretion of diffuse matter and minor clumps from the 
perturbation wings (details and further references are 
given in Appendix Al). The two stages are separated 
by the redshift Zt ~ 0.5 — 1 when the circular veloc- 
ity at the virial boundary attains its maximal value; 
this epoch leaves a clear imprint in the halo 'concentra- 
tion' parameter c = i?/r_2 (the ratio of the virial radius 
to the reference radius in the halo bulk where the DM 
density slope equals —2), that grows after zt following 
c(z)«3.5(l4-zt)/(l + z). 

After Zt the halos attain a quasi-static equilibrium de- 
scribed by the Jeans equation; the explicit solutions ('a- 
profiles', with a — 1.27 in rich clusters) for the den- 
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sity p{r) and the gravitational potential $(r) are given 
by Lapi & Cavaliere (2009a, b) and recalled in Eq. (B3). 
Note that the physical scales including r_2 are modu- 
lated by c. 

We will fulfill our purpose with the use of two main 
tools: basic entropy patterns, and the entropy-based 
equilibrium condition expressed by our Supermodel. 

2. ENTROPY PATTERNS 

The basic ICP entropy run we expect to apply through- 
out the ICP can be rendered as a central level fee connect- 
ing to a rising ramp with slope a toward the outer value 
k2 in the form (see Tozzi & Norman 2001, Voit 2005) 

k{r) = kc + k2 {r/Rf , (2) 

illustrated in Fig. 1 (red solid line). Next we discuss the 
physical origin of such a minimally structured distribu- 
tion. 

2.1. Central entropy 

In the central range r < 2 x 10^ kpc the entropy is ini- 
tially set at levels fee ^ 10 keV cm^, not much exceeding 
the levels fci prevailing in the IGM (see Ryu et al. 2008, 
Nicastro et al. 2010). This is because during the initial 
fast collapse the temperatures in the virialized core are 
high, at fcsT « GmpM{< r)/10r ^ a few keVs, but 
the ICP is dense at n ~ 10"'^ cm~^, in step with the 
general overdensities S p/p > 2 x 10^ over the average 
environment. 

Such entropy levels are eroded or even erased away 
following ds/dt = —s/tc, due to the radiative cool- 
ing by bremsstrahlung (increasingly dominating over 
line emission for fc^T > 2 keV) that makes up the 
observed X-ray emissions; the associated timescale for 
a single-phase ICP (see Sarazin 1988) reads tc ~ 
30 {ksT /kcV)^/^ {n/10-^ cm"-"')-! Gyr. Thus cooling 
may be slow and little relevant in the low-density out- 
skirts, but is speeded up in the dense central ICP, so 
that within some 5 Gyr the levels kc are depressed from 
10^ keV cm^ down to ~ 10^ keV cm^. Wherefrom 
cooling becomes so fast as to match the dynamical times 
~ 10^^ Gyr, to the effect of impairing the thermal pres- 
sure support; the process is even faster in multi-phase 
ICP with a considerable cold component. 

This leads to ICP condensation and to cooling faster 
yet, so starting up an accelerated settling to the cluster 
center and onto the central galaxies (the classic 'cooling 
catastrophe'), were it not for renewed energy injections 
(see Binney & Tabor 1995, Cavaliere et al. 2002, Voit & 
Donahue 2005, Tucker et al. 2007, Hudson et al. 2010). 
These occur when the accretion reaches down into the 
galactic nuclei and onto their central supermassive black 
holes, to trigger or rekindle a loop of intermittent star- 
bursts and AGN activity; in the form of gentle bubbling 
or moderate outbursts over some 10^^ Gyr, this can sta- 
bilize the time-integrated kc at levels of ^ 10^ keV cm^. 
Such an enticing scenario is discussed, among others, by 
Ciotti & Ostriker (2007), McNamara & Nulsen (2007), 
Conroy & Ostriker (2008), Churazov (2010). In sum, a 
cool core constitutes an attractor for the thermal state 
of the central ICP. 

On the other hand, kc may be raised up to levels 
of several 10^ keV cm^ when substantial energy injec- 
tions AE occur into the ICP from violent outbursts of 



AGNs in central galaxies or even more from deep merg- 
ers. These injections drive through the central ICP a 
blastwave bounded by a leading shock with Mach num- 
ber gauged by the relation 1 + AE/E in terms 
of the ICP thermal energy E ^ 2x 10^^ (fcijT/keV)^/^ 
erg (see Lapi et al. 2005, their Fig. 7); a strong shock 
with Ai'^ > 3 would require injections AE >2E, i.e., a 
few tens of keVs/particle. This may be the case for deep 
major mergers, while it is hardly matched by an AGN 
powered by a supermassive black hole up to 5 x 10^ Mq 
with only some 5% of the energy discharged effectively 
coupled to the ICP (see Lapi et al. 2005). 

Blasts that preserve the overall equilibrium may still 
leave a long-lasting imprint onto the central ICP in 
the form of an entropy addition spread out to a radius 
r/ « 10^ kpc where the blast has expanded, stalled and 
degraded into sound waves. A handy representation (see 
Fusco-Femiano et al. 2009; see also Appendix B) of the 
entropy distribution in such conditions is still given by 
Eq. (2) for r > r/, while a roughly constant level 

k{r) = kc (3) 

applies for r < r/, as illustrated in Fig. 1 (orange dashed 
line). We shall see that conspicuous central wiggles may 
then appear in the radial temperature profiles if these 
imprints survive spherical averaging; such features will 
persist over timescales longer than the blast transit time 
^ 0.3 Gyr if shorter than the cooling time ~ 5 Gyr. 

Stronger if rarer energy injections with AE ^ E can 
be produced as head-on major mergers following the halo 
bulk collapse (see McCarthy et al. 2007, Norman 2011) 
deposit at the center large energies of several tens keV per 
particle, and entropy levels » 10^ keV cm^; these trigger 
conditions of severe disequilibrium such as in A754 (see 
Macario et al. 2011) and A2146 (see Russell et al. 2010), 
or outright disruption like in MACS J0025.4-1222 (see 
Bradac et al. 2008), or 1E0657-56 (the 'Bullet Cluster', 
see Clowe et al. 2006). 

2.2. Outer entropy 

Supersonic inflows of external IGM drive at the cluster 
boundary R ~ Mpc strong shocks intertwined into a web 
or layer located at i?^ « i? where accretion feeds on 
filaments (see Lapi et al. 2005, Voit 2005). These shocks 
are effective in thermalizing a considerable fraction of the 
specific energy vf that the IGM gains as it infalls from 
the 'turnaround' radius i?ta to the virial R « Rta/2, 
to the effect of providing substantial temperature jumps 
T2/T1 from the IGM values. These jumps grow with 
the Mach number squared Al^, whilst the density jumps 
n^jnx saturate to 4, and the post-shock kinetic energy 
v\lv\ decreases with A4~^ in the shock restframe, cf. 
Appendix A and Fig. A2. As a result, soon after the 
cluster formation large thermal energies are deposited 
in the thin ICP at the boundary, with densities close 
still to the IGM's; there the entropy levels reach up to 
5 X 10^ keV cm^. 

The bearing of these issues to the ICP physics is fo- 
cused from the expression derived by Cavaliere et al. 
(2009) for the value of the entropy slope an at the bound- 
ary 

ail = 2.5- 0.5 6fl . (4) 
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This value (clearly smaller than 2.5) sensitively depends 
on the ratio bn = /imp v^/kBT2 > 1 of the potential to 
the ICP thermal energy at r = i? (see Lapi et al. 2005). 
Values an « 1.1 are obtained soon after the bulk collapse, 
when the inflow is still sustained and strong shocks fully 
thcrmalizc the infall energy = 2 A"I> into three degrees 
of freedom, and produce postshock temperatures kBT2 ~ 
HmpVi/3 (for closer evaluations see Appendix A3-A4). 
On expressing the potential drop from the turnaround 
to the shock in the form A^/vj^ « 0.57 (sec Cavaliere 
et al. 2009), the standard values « 3vjj/2A^ w 2.7 
and a w 1.1 are obtained (Tozzi & Norman 2001). 

Eq. (4) is derived as the current boundary value ur 
for a, but it clearly yields also the running slope a(r) in 
the middle range on considering that — in the absence 
there of energy sources — the entropy will be conserved 
and stratified at the values previously produced when 
the boundary was just at r. In other words, the radial 
entropy distribution preserves the memory of the past- 
time development. 

As the cluster outskirts grow farther out, the inflows 
slow down considerably, and do so especially at low z; 
this straightforwardly is to occur when the accretion is 
drawn from the tapering wings of a DM perturbation 
over a background lowering under the accelerated cosmic 
expansion. Thus the potential drop A$ becomes shal- 
lower (see Appendix Al; also Lapi et al. 2010) while the 
shocks outgrow R, to the eff'ect of weakening the shock 
jumps and lowering T2 toward the external value; Ti. As 
a result, bfj grows and a decreases toward zero. 

A handy representation (see Lapi et al. 2010; see also 
Appendix B) of the ensuing entropy distribution is still 
given by Eq. (2) inside r < with rt R/3 (to be 
discussed in § 4), while 

k{r) = k2 (r/i?)«+"' e"' (^-^y^^" (5) 

applies for r > r;,, as illustrated in Fig. 1 (blue dotted 
line). This expression describes a simple linear decline of 
the slope a(r) with a gradient a' = (a — an)/{R/ri, — 1) 
from the inner value a ^ 1.1 to the outer value aji < a. 

Such an entropy bending takes place on the timescale 
set by the outskirts development, when the DM halo 
grows its concentration to values c > 6 from the initial 
values c « 3.5 set soon after the bulk collapse at zf, e.g., 
for a cluster collapsed at Zt ^ 1 and observed at 2; « 0.15 
the time elapsed amounts to 6 Gyrs. 

In sum, the outer ramp flattens and bends over a 
timescale of several Gyrs, while the central level kc is 
eroded away by radiative cooling. These two changes 
are independently driven at far apart locations by quite 
different processes; what they have in common, though, 
is their progress in time. So one expectation from our 
picture is that they should take place together as the 
structures age, a main feature in our cluster classiflca- 
tion of § 5. 

3. THE ENTROPY-BASED EQUILIBRIUM 

The entropy-based equilibrium of the ICP within the 
DM gravitational wells is constituted by our Supermodel 
(see Cavaliere et al. 2009 and Fusco-Femiano et al. 
2009), with the related straight algebra recapped in Ap- 
pendix B. 



3.1. Thermal support 

There we recall that the linked radial proflles of tem- 
perature and density read 

kBT{r) = n2/3(r) k{r) cx k^/^r) [1 + 2/5 bnl{r)] , (6) 

having used the shorthand T(r) = 
J^^j^dx[vl{x)/vj^][k{x)/k2]~^^^ /x in terms of the 

circular velocity v"^ = GM/R (see also Appendix B, 
below Eq. B3). As discussed by Cavaliere et al. (2009), 
in the outskirts I is small compared to 1 and the whole 
factor in square parenthesis behaves like (r/i?)"^*"/^; 

on the other hand, at the center the integral T oc fee 
dominates over 1 and scales inversely with the central 
entropy level kc- 

These temperature and density proflles provide the vol- 
ume emissivity for bremsstrahlung, proportional to 

ocn2(^)Ti/2(r) oc fc-9/i°(r) [1 + 2/5 6^1(0]^/' ; 

this constitutes the basis for computing (after spectral- 
bandpass windowing and projection) the X-ray observ- 
ables, namely, the surface brightness and the emission- 
weighted temperature; full expressions are given in Ap- 
pendix B. 

We stress that all these profiles for n(r), T(r), and 
Sx (r) are linked together by the underlying entropy dis- 
tribution. For a relevant example, Eq. (6) yields the 

central scaling laws Tc cx fc°-^^ and k^ , see 

Cavaliere et al. (2009); thus when kc is low the temper- 
ature will dip and the associated emissivity will rise to- 
ward the center, features that constitute the marks of the 
conventional cool-core designation. On the other hand, 
high kc produce flat emissivity proflles together with a 
wide temperature plateau, typical of the conventional 
non-cool-core designation. Moreover, the central cool- 
ing time in single-phase equilibrium may be expressed 
in terms of the entropy level kc only, to read simply 
tc » 0.5 (fcc/15keV cm^)^-^ Gyr; this implies that high 
levels of kc ~ 150 keV cm^ require timescale of order 8 
Gyr to be eroded. In the outskirts, instead, the scaling 
T(r) oc r^l^o-R-"^ holds, showing that when k(r) is bent 
down with a/j <C 1 the temperature will fall steeply out- 
wards; in simple terms, the profile T(r) cx k(r)r?l'^(r) 
will follow n^/'^(r) or steeper when k(j) is nearly con- 
stant or even bent down. Meanwhile, the brightness will 
be flatter at intermediate radii (see Fig. 2), and con- 
stitutes a simple pointer toward interesting temperature 
and entropy distributions (see Cavaliere et al. 2011). 

An observable independent of X-ray data is provided 
by the SZ effect (Sunyaev & Zel'dovich 1972); the radial 
profile of its strength parameter is proportional to the 
thermal ICP pressure and writes as 

y(r) cx n(r) r(r)cx [1 + 2/5 &flX(r)]'/' . (8) 

We stress that the Supermodel implies a nearly universal 
pressure profile (and correspondingly for the SZ effect), 
since the entropy radial dependence is encased into the 
slowly varying factor X(r); this is the ultimate origin for 
the approximate invariance of the pressure profile derived 
from the X-ray data by Arnaud et al. (2010). Using the 
inner scaling of n and T with kc we find that the scaling 
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y cx holds, implying that higher values of y cor- 

respond to lower kc- At the other end, in the outskirts 
y oc r^""^'"^ applies, implying sharper declines in clus- 
ters with shallower entropy ramps. Thus the (projected) 
SZ effect provides a direct probe of the entropy levels 
throughout a cluster, and so an independent way for clas- 
sifying HE and LE types from Plnnck (see Aghanim et 
al. 2011) and from ground-based instruments. 

3.2. Turbulent support 

As argued above, the conditions of low entropy pro- 
duction are related to mildly supersonic inflows and weak 

boundary shocks with decreasing Mach number Ai^ < 3; 
we stress that in turn they are conducive to trigger outer 
subsonic turbulence developing under the drive of rela- 
tively more inflow energy V2 / vf oi A4~'^ seeping throiigh 
the weaker shocks (see Cavaliere et al. 2011; also Ap- 
pendix A4 for details). The turbulent contribution to 
equilibrium is conveniently described in terms of the 
ratio S = Pnth/pth of the turbulent to thermal pres- 
sure. The boundary normalization is consistently set by 
Sr oc V2/vf, while the shape S{r)/SR of its inward decline 
on a scale ^ ~ 10^ kpc is provided by the classic cascade 
from large 'eddies' at the macroscopic coherence length, 
fragmenting to small eddies where dissipation becomes 
effective (see Kolmogorov 1941, Monin & Yaglom 1965; 
see § B3 in Appendix B for details). 

In fact, it turns out that the total pressure pth +Pnth = 
Pth (1+^) can be straightforwardly included in the hydro- 
static equilibrium solved by the Supermodel; the result 
can be described simply in terms of Eq. (6), with T and 
k replaced everywhere (including I) by T = T {1 + 6) 
and by 

~k = k{l + S). (9) 

The underlying rationale is that turbulent eddies add to 
the microscopic thermal degrees of freedom in dispersing 
and ultimately dissipating the inflow kinetic energy t;| 
seeped through the shock. 

While turbulence is stirred, the thermal pressure re- 
quired for overall support in the given DM gravitational 
potential well is decreased. If turbulence were not ac- 
counted for, the overall masses estimated from X rays 
would tend to be negatively biased compared to the grav- 
itational lensing measurements (Nagai et al. 2007, Lau 
et al. 2009, Meneghetti et al. 2010, Kawaharada et al. 

2010) . Meanwhile, the intensity parameter of the volume 
thermal SZ effect y{r) is lowered relative to the pure ther- 
mal equilibrium expression Eq. (8) by an explicit factor 
l/{l + 6), adding to small corrections to the integrand in- 
side X. Note that such a straightforward lowering is con- 
siderably stronger than may result from any reasonable 
ion-electron disequilibrium at the shock (see the accurate 
estimates by Wong & Sarazin 2009). Thus SZ effect can 
also provide a direct probe of a low thermal pressure, 
which implies a considerable turbulent component in the 
cluster outskirts for attaining equilibrium (see Cavaliere 
et al. 2011). The dearth of outer thermal pressure is 
indicated by stacked WMAP data (see Komatsu et al. 

2011) ; the contribution to such conditions from LEs and 
HEs will be discussed in a forthcoming paper. 



4. TOWARD A CLUSTER LIBRARY 

We aim at constructing flrst a library of clusters from 

extended circularly-averaged data; on this basis we aim 
at introducing a physically meaningful cluster classifica- 
tion scheme and at discussing the connection of the ICP 
thermal state to the DM halo development. This requires 
robust fits to the X-ray observables from linked, consis- 
tent profiles of density and temperature, as to pinpoint 
the few independent parameters governing the ICP en- 
tropy distribution. Speciflcally, we adopt the following 
strategy. 

The entropy-based picture of § 2 (illustrated in Fig. 1) 
suggests the basic entropy distribution of Eq. (2), consti- 
tuted by a central level kc going into a ramp rising with 
slope a « 1 toward the outer value ^2 . Two relevant and 
alternative variants may apply: the central floor kc ex- 
tends out to a radius r j and is angled there to the outer 
ramp, see Eq. (3); beyond a radius Vb the ramp bends 
over to a shallow slope ai? <C 1 joining a low boundary 
value k2, see Eq. (5). The basic distribution provides a 
baseline with a minimal number of parameters; the flrst 
variant is convenient when central temperatures are high, 
but wiggles stand out, while the second applies to cases 
with low central temperature and steep temperature de- 
cline into the outer region. 

We first insert the basic entropy distribution Eq. (2) 
with free shape parameters kc and qr into the Su- 
permodel Eq. (6), and derive the radial profiles of 
density and temperature; the Supermodcil algorithm is 
made available at the URL http://people.sissa.it/ 
~lapi/Supennodel/. We then perform a fit to the pro- 
jected, emission-weighted temperature data, using a mul- 
tiparametric minimization procedure (e.g., MPFIT by 
Markwardt 2009), and derive the temperature scale T2 
at the boimdary (sec Eq. B7). Finally, we fit the sur- 
face brightness including the bandpass correction (that 
requires T2), and derive the scale n2 at the boundary 
(see Eq. B8). Thus we obtain also the entropy normal- 
ization k2 = kBT2ln2^ at the boundary, to complete the 
entropy distribution. 

When the value of a fit turns out to be large, we 
proceed to insert in the Supermodel the variant entropy 
distributions given by either Eq. (3) or Eq. (5); this adds 
a further parameter, either ry- (for A644 and A2256) or 
rb (for A1795, PKS0745, A2204, A1413). The statistical 
signiflcance of the added parameter is corroborated by 
relevant improvements in the reduced values (see Ta- 
ble 1). This is further checked with the F— test, which 
yields a significance level larger than 98% but for A1413 
where it is larger than 96%. 

We stress that our Supermodel fits are performed over 
the whole radial range covered by the current X-ray data. 
In a number of clusters observed by Suzaku (e.g., A1795), 
the X-ray data extend out to approach the virial radius 
R; in other instances observed by XMM-Newton (e.g., 
A1656) the data are more limited (around -R500), and 
the outer parameters are provided by the Supermodel 
upon extrapolation, implying larger uncertainties. In the 
case of PKS0745 the uncertainty is particularly large due 
to discrepancies between the XMM-Newton and Suzaku 
datasets, see the discussion by Eckert et al. (2011). 

We note that the shape parameters kc and ur may 
be determined from fitting either the temperature or the 
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brightness profile; the results arc consistent within the 
respective uncertainties, but the value derived from the 
former is to be preferred whenever extended, high-quality 
data are available (with the due caveats discussed by 
Eckert et al. 2011 as to anomalous background, and by 
Simionescu et al. 2011 as to effects of possible dumpiness 
in one sector of the Perseus cluster), since the tempera- 
ture dominates the brightness in the entropy expression 
k oc T^/<5 5.-1/3 (-ggg Cavaliere et al. 2005). 

On the other hand, fits to the X-ray brightness can 
also provide the DM concentration c = Rjr^i that en- 
ters the Supermodel formalism through (r) , while the 
outer scale R is provided by independent observations 
such as the red-sequence termination or gravitational 
lensing. We note that the determination of c is mainly 
based on the outer brightness data, so is closely indepen- 
dent of the inner entropy distribution. We stress that 
the Supermodel leads to a fast yet robust evaluation of 
c from X rays only, with results less biased than gravita- 
tional lensing by prolateness effects (discussed, e.g., by 
Corless et al. 2009). 

We also note that the few parameters entering the en- 
tropy distribution are calibrated from fitting with the 
Supermodel the observables directly expressed in terms 
of the radial profiles of n(r) and T'(r), with no need for 
delicate data deprojections (discussed by Yoshikawa & 
Suto 1999, Cavaliere et al. 2005, Croston et al. 2006, 
Urban et al. 2011). 

The ICP parameters so derived are used here to build 
up the library of 12 clusters presented in Table 1. Ten of 
these have been analyzed by us in previous works (Cav- 
aliere et al. 2009, Fusco-Femiano et al. 2009, Lapi et al. 
2010), while here we add A399 and refine the analysis of 
A2218. Examples of Supermodel fits are illustrated in 
Fig. 2. 

One may ask to what extent the entropy distributions 
derived from the Supermodel might depend on the un- 
derlying assumptions. These concerns are swept away by 
Fig. 3 where we compare, in the radial range r > 0.2 i? 
where cooling is negligible, the entropy distributions de- 
rived from our Supermodel analysis of the 12 clusters 
listed in Table 1, with the outcomes of the nonradiative 
hydro-simulations by Burns et al. (2010) for a sample of 
24 relaxed massive clusters. Our results are seen to be 
consistent with the simulation outcomes and their vari- 
ance, that grows wider into the outskirts. 

Moreover, such comparison implies that throughout 
most of the cluster volume the Supermodel results are 
robust against the assumptions of spherical symmetry, 
hydrostatic equilibrium, and purely smooth accretion. 
In fact, in the inner regions merger-related geometrical 
asymmetries are smoothed out on a crossing timescale, 
shorter than the time required by cooling to erase en- 
tropy excesses of ~ 10^ keV cm^. In the middle regions, 
approximately spherical symmetry of the ICP is indi- 
cated by various simulations (e.g., Lau et al. 2011). In 
the outer regions, the accretion is dominated by minor 
mergers or truly diffuse matter, as shown in detail by the 
simulations of Wang et al. (2011, see their Fig. 7). All 
that explains why the snapshots provided by our Super- 
model fits to the X-ray data are so robust. 



5. A CLUSTER GRAND DESIGN 

We have parted the 12 clusters listed in Table 1 into 

two main blocks on the basis of their k,, values being of 
the order of a few 10^ or a few 10^ keV cm^; within each 
block, we have ordered the chisters on the basis of their 
values. It is easily perceived that the two main blocks 
are also parted as to the values of their DM concentration 
c. The ordering indicates correlations between these ba- 
sic physical parameters quantified in § 5.1, and suggests 
an evolutionary trend linking the ICP thermal state with 
the DM development to be discussed in § 5.2. 

5.1. Correlations 

In the top panel of Fig. 4 we illustrate the central en- 
tropy level kc vs. the entropy slope a in the cluster bulk 
at i?5oo « 0.5 R. For cool-core clusters our results (blue 
dots) from the Supermodel analysis compare well as to 
central values and their uncertainties with the sample of 
relaxed, mostly cool-core clusters by Pratt et al. (2010; 
green squares). It is seen that as to the average data 
values, kc correlates poorly with a at i?5oo, as quantified 
by the low value of the Spearman's rank correlation co- 
efficient p « 0.27 (cf. Lupton 1993) both for our and the 
above authors' samples. 

In the middle panel of Fig. 4 we illustrate the cen- 
tral levels kc vs. the outer slopes an. As to the latter we 
find values close to a for non-cool-core clusters (red dots), 
while appreciably lower for cool-core clusters (blue dots). 
It is seen that an correlates well with kc as to the average 
values; this is quantified by the value of the Spearman's 
rank correlation coefficient p » 0.64. On the other hand, 
the often large uncertainties in or (related to true un- 
certainties in the outer X-ray data) and especially in kc 
(related also to inner physical complexities) will blur the 
correlation. 

We test to what degree this occurs by running 10^ 
Monte Carlo simulations, randomly sampling values of 
kc and an from Gaussian distributions around their av- 
erages, with widths given by their formal 1 — a uncertain- 
ties in both variables; with this conservative treatment 
the average Spearman's coefficient is lowered to p « 0.44, 
corresponding to a 9% probability for chance occurrence 
of the correlation. In addition, we compute that on sta- 
tistical grounds the probability of 'outliers' (objects with 
kc > 30 keV cm^ and an < 0.6) is 5% on average, with 
a formal standard deviation of 22%; this implies that on 
doubling the size of present sample to 24 objects, one 
should expect from 1 to 7 outliers. The above outcomes 
motivate us to investigate in § 5.2 whether a physical 
basis underlies the apparent dearth of clusters with high 
central entropy levels kc > 10^ keV cm^ and low outer 
entropy production an < ^■ 

In the bottom panel of Fig. 4 we illustrate the outer 
slopes an vs. the concentration parameter c, derived 
with the Supermodel. We find that low values of an 
correspond to high values of c, which mark a long life- 
time from the formation zt to the observation redshift 
^; « following c « 3.5 (1 -I- ^t)/(l + z). Such an anti- 
correlation between an and c is highly significant as for 
the average data values, with a Spearman's coefficient 
p « —0.79; on the other hand, the additional uncertain- 
ties that also affect c lower it to the conservative value 
p w —0.46. Again, this outcome stimulates us to inves- 
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tigate any physical dearth of ckistors with high concen- 
tration c > 6 and steep outer slopes aj^ > 1, and offers a 
pattern to be confronted with future real and/or virtual 
datasets. 

5.2. Classes 

Guided by the above discussion, we submit that all 
our clusters may be parted into two main classes, defined 
on the basis of low entropy (LE) or high entropy (HE) 
prevailing not only in the inner region but also throughout 
the ICP. 

• HE clusters, featuring high entropy throughout the 
ICP; that is, featuring not only a central level kc ss 
3 X 10^ keV cm^, but also a vc;ry high boundary level 
A;2 ~ 3 — 5 X 10'^ keV cm^ corresponding to a steep en- 
tropy ramp with a > 1 throughout the outskirts. The 
high values of kc yield a monotonic temperature pro- 
file T(r) throughout, slowly declining from the central 
plateau into the outskirts. We stress that our class defi- 
nition includes not only a central non-cool-core state as 
in the designation introduced by Molendi & Pizzolato 
(2001) and pursued by Leccardi et al. (2010), but also 
an associated high level of outer entropy production. We 
propose that the association arises due to the young age 
of the containing DM halos, marked by low values of the 
concentrations c « 4, with a lifetime too short for central 
entropy to be erased away and any entropy bending to 
be effective in the outskirts. 

• LE clusters, featuring low entropy throughout the 
ICP; this includes both a low central baseline kc < 30 
keV cm^ and a moderate outer level k2 < lO'' keV cm^, 
so as to imply a ramp bending outwards of r^/i? > 0.3 
toward a{r) < 1 (see also Hoshino et al. 2010); the out- 
come is a low central value of T and a peak of T(r) at 

^ 0-2 followed by a steep decline outwards, partic- 
ularly effective at low z (e.g., A1795). Om class defini- 
tion includes not only a central cool-core state as in the 
standard designation, but also an associated low level of 
outer entropy production. We propose that the associa- 
tion low kc — shallow an is to be traced back to the long 
lifetime of the containing DM halos, marked by high val- 
ues of the concentrations c > 6. We relate such a late 
stage in the outskirts development to dwindling inflows 
that cause weaker boundary shocks with < 3, low 
entropy production and a substantial fraction of kinetic 
energy left over to drive outer turbulent eddies. 

The low kc levels proper to LEs are driven by cooling 
timescales tc shorter than the halo dynamical age marked 
by c. In fact, the divide between LEs and HEs is around 
kc ~ 150 keV cm^ corresponding to a cooling time tc ^ 8 
Gyr (e.g., the lapse between 2 « 1 and z « 0.1); after 
this, fast cooling leads to an accelerated progress toward 
kc levels lower yet. Eventually, however, the levels of 
kc are likely to be stabilized by two additional physical 
processes, i.e., intermittent AGN activity and impacts 
of deep major mergers; two modes are suggested by the 
broad, possibly double-peaked distribution for the num- 
ber of clusters with given kc, as observed by Cavagnolo 
et al. (2009) and Pratt et al. (2010), and discussed by 
Cavaliere et al. (2009). 

The relationship between the classes is illustrated in 
the evolutionary chart of Fig. 6, that represents our clus- 
ter Grand Design. This envisages clusters mainly born 
in a HE state of high entropy, dominated by the fast vi- 



olent collapse of the halo bulk and related strong shocks 
in the infalling gas. Subsequently, on a timescale of sev- 
eral Gyrs they progress toward a LE state since both the 
central entropy is lowered by radiative cooling, and the 
outer entropy bends over because of the weakened shocks 
and tapering entropy production. The Grand Design en- 
visages that in a number of cases such a sequence may be 
halted within a few Gyrs and reversed by late, trailing 
deep mergers which remold any nascent cool-core and 
rejuvenate the central ICP into a higher entropy state. 

In fact, these clusters with ICP lingering in such an 
intermediate state may be conveniently ranked in a sub- 
class labeled HE, marked out from basic HE by a wiggled 
central temperature profile. We have shown (see discus- 
sion by Fusco-Femiano et al. 2009) that such profiles 
obtain whenever the central entropy features a floor kc 
extended out to a radius rf ^ 10^ kpc; correspondingly, 
the central brightness features a particularly flat proflle. 
We recall from § 2 that such entropy additions are likely 
imprinted by a blastwave with Mach numbers A^^ > 3 
launched outwards by a head-on impact of a deep merger. 
When the blast has stalled around rf and the overall 
equilibrium in the ICP is recovered, the central entropy 
is still enhanced to levels up to kc ^ 10^ keV cm^, and 
so is immune to subsequent, weaker AGN-driven blasts. 

Such a HE morphology turns out to occur not only in the 
two cases listed in our Table 1, but also in several more 
instances of the kind illustrated by Rossetti & Molendi 
(2010), close to 50% of their non-cool-core clusters; thus 

we propose the HEs to deserve a subclass status. 

Our interpretation of the HE morphology relates the 
size Tf to the epoch of the merger responsible for the en- 
tropy input; such an epoch is expected to be in between 
the blast transit time rf/Mcg ~ some 10~^ Gyr and 
the several Gyrs taken by radiation to erode the floor, 
or by central turbulence to blur it. Such a timing en- 
sures an accurate description of the ICP thermodynam- 
ics by the Supermodel based on hydrostatic equilibrium. 
Note also that the ICP attains its equilibrium somewhat 
faster than the DM does (see Ricker & Sarazin 2001, 
Lapi et al. 2005), while the circularized data (see Snow- 
den et al. 2008) tend to smooth out limited deviations 
from spherical hydrostatics (for a detailed discussion see 
Fusco-Femiano et al. 2009). 

6. PREDICTIONS FROM THE GRAND DESIGN 

Here we present other, specific predictions derived from 
our cluster Grand Design. 

• We expect the HE clusters to feature outer profiles, 
r(r) declining mildly to a boundary value T2 still sus- 
tained. We illustrate in Fig. 5 our prediction for such a 
mild decline in the two HE clusters A1656 and A2256. 
compared with the ciirrently limited data. Note that 
for A1656 our Supermodel fit has only used the data by 
Snowden et al. (2008) out to r « i?/3; our outer pre- 
diction agrees with the recent data by Wik et al. (2009) 
extending out to R/2. 

• We expect LE clusters to feature at low z particu- 
larly small values of kc and sharply bent outer entropy 
profiles. The latter yield steeply declining T{r) profiles, 
as supported by the Suzaku observations of a few clus- 
ters like A1795; a similar recent case may be constituted 
by A2142, see Akamatsu et al. (2011). Low SZ signals 
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and increasing contribution of outer turbulent support 
are also expected; relatedly, in these clusters the mass 
reconstructed from X-ray observations will show system- 
atic deficits relative to the gravitational lensing result 
(see Cavaliere et al. 2011). 

• We expect at higher z a lower fraction of LEs, re- 
flecting the main evolutionary trend from HEs to LEs 
envisaged by our Grand Design; this is consistent with 
the evidence by Santos et al. (2010) based on observ- 
ing the average surface brightness up to redshift 2: « 1.3. 
When observations of very low surface brightness will 
become feasible, we expect steeper brightness profiles 
and a milder temperature decline in the outskirts to 
loom out (see § 3.1), as for such high-z LEs the cos- 
mology/cosmogony had not time enough to sharpen the 
outer entropy bending. 

Our picture envisaging low or high entropy levels to 
hold throughout the ICP is consistent with the present 
dearth of the following pairings: nearby clusters with 
low kc levels and high an values (that would be located 
in the upper left strip of the kc — an plane in Fig. 4); 
clusters with high kc > 10^ keV cm^ and low < 1 
(that would be located in the lower right corner of the 
kc — o-R plane in Fig. 4); inner temperature wiggles in 
highly concentrated chisters with c > 6. Wider libraries 
based on extended, high-quality temperature data will 
allow testing the above predictions. We add that the 
Supermodel predicts the projected SZ effects (otherwise 
closely universal, § 3.1) to differ from the HE to the LE 
cluster population, with the latter featuring steeper pro- 
files in the outskirts; in fact, as stated in § 3.2 in LEs 
outer turbulence is expected to contribute substantially 
to the equilibrium, so as to lower by 1 + 15 the thermal 
SZ effect. 

7. DISCUSSION AND CONCLUSIONS 

We have seen how the cluster thermodynamical state of 
nearly relaxed clusters can be probed by means of linked, 
robust profiles of density and temperature in the ICP 
derived from extended X-ray data. We have carried out 
such a task on using the Supermodel formalism with its 
few, intrinsic parameters that modiilate the underlying 
distribution of the specific entropy. On this basis, we 
have grouped the rich clusters analyzed here into two 
main classes, on account of low (LE) or high entropy 
(HE) prevailing throughout the ICP. 

Such classes constitute thermal conditions with long 
persistence. In fact, HEs with their hot atmospheres 
are stabilized by long cooling times and stubborn resis- 
tance to supersonic flows; LEs are likely stabilized by in- 
ner AGN energy injections, while their outskirts evolve 
slowly as the inflows across the boundary decrease to- 
ward low z. The main overall evolutionary course pro- 
ceeds from HEs to LEs due to erosion of central entropy 
by cooling, and to reduced production of outer entropy 
by weakened accretion shocks. 

However, such a course may be interrupted or even 
reversed by large entropy injections from major merg- 
ers, particularly frequent at high z; thus an HE thermal 
state sets in, marked out by wiggles in the central radial 
temperature distribution. These are interpreted in the 
Supermodel framework in terms of a sharp entropy floor 
extending out to r/ > 50 kpc and with levels kc around 



200 keV cm^. We consider these as intermediate objects, 
constituting a subclass contiguous to, and blending into 
the HE main class. 

Our overall picture derived from the snapshots con- 
densed in Table 1 relates the ICP thermal state to the 
DM halo development stage, in the form of an inverse 
correlation between the outer entropy slopes a r and the 
halo concentrations c. While the outer ICP thermody- 
namical age is signaled by the former, the DM dynam- 
ical age is marked by the latter, specifically in terms of 
c(z) « 3.5 (1 -I- Z()/(l -I- z) increasing from the formation 
zt to the observation redshift z. 

We interpret such a correlation as follows. The LEs 
are associated with high-c halos, old enough as to al- 
low the ICP to be affected by deep radiative erosion of 
their central entropy (producing low kc values) and by 
reduced entropy production in the outskirts (shallow ur 
or low /C2); the latter effect inescapably depends on large- 
scale cosmogonical/cosmological evolution, and at given 
c is more pronounced at low z when reduced accretion 

is most effective. Conversely, the HEs (and HEs) are 
associated with young halos of low c < 4. We stress 
that such central and outer ICP thermal evolutions are 
independently driven at far apart locations by different 
processes; what they have in common, though, is their 
largely parallel progress over comparable timescales of 
several Gyrs. 

On the other hand, a reasonable amount of variance in 
central entropy and in outer bending may produce some 
intermediate instances; one such case is constituted by 
A1689 at z w 0.18, with its still rather high kc level 
and intermediate values of its outer entropy slope. As a 
matter of fact, considerable variance around the average 
picture will be caused by the well-known scatter in the 
birth and development of cosmic structures; this affects 
both the halo collapse redshifts Zj (e.g.. Bullock et al. 
2001, Wechsler et al. 2006, Klypin et al. 2011), and 
the subsequent merging histories (e.g., McCarthy et al. 
2007. Fakhouri et al. 2010). In particular, A1689 with its 
mass Af s::^ 1.3 X 10^"'' Mq constitutes a well-studied case 
of high zt ~ 2.5 as inferred from its high concentration 
c > 10 (see Appendix A; also Broadhurst et al. 2008, 
Lapi & Cavaliere 2009b). 

Another source of variance is constituted by the cluster 
environment; in particular, adjoining filaments with con- 
trasts 5p/p ^ 5 will enhance diffuse accretion in a rich 
ambient like a supercluster, so as to delay weakening of 
shocks and onset of turbulence. This may be the case 
with one sector out of four in A1689 (see Kawaharada et 
al. 2010, Molnar et al. 2010) and with A2199 (see Rines 
et al. 2002), implying the spherically averaged values of 
a to be higher than in standard LE clusters; the opposite 
holds true for cluster sectors facing voids. 

Additional variance might arise by cold subclumps in 
sectors of the nearby Virgo and Perseus clusters; this 
would bias high the surface brightness and the apparent 
baryonic fraction (see Ettori et al. 1998 vs. Simionescu 
et al. 2011, Urban et al. 2011). On the other hand, 
such features do not appear to affect most of the clus- 
ters collected in Table 1, including instances with steep 
temperature decline and flat entropy distributions like 
A1795; as for the latter, our Supermodel yields an outer 
baryonic fraction bounded by 0.14. Similar values have 
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been recently inferred from aimed Suzaku observations 
of A2142 (see Akamatsu et al. 2011). 

In LE clusters we expect outer turbulence related to 
compressive modes to develop under the drive of kinetic 
energy increasingly seeped through weakening shocks 
(Cavaliere et al. 2011). Inner turbulence, on the other 
hand, is likely stirred in HE clusters by shear motions 
associated with the mergers' wakes (e.g., lapichino et al. 
2011). These motions are widely held to accelerate elec- 
trons in situ up to Lorentz factors 7 ^ 10**; the electrons 
energize strong radiohalos by their synchrotron radiation 
in cluster-wide magnetic fields of a few /iGs, with electron 
lifetimes under 1 Gyr (see Ferrari et al. 2008, Feretti et 
al. 2011, Brunetti 2011), shorter than the thermal cool- 
ing times around 5 Gyr for the center of an HE cluster. 
On the other hand, in the process of cooling toward an 
LE state with levels kc < 50 keV cm^ the core becomes 
sufficiently cold as to be sensitive even to lesser merg- 
ers. Then temperature wiggles and radiohalos may form 
together, but the latter will fade much sooner than the 
former can be eroded away (see Buote 2002, Brunetti et 
al. 2009; also Rossetti et al. 2011); as a result, we expect 

more HE clusters than radiohalos. 

In this paper, we have shown how entropy offers a key 
to detailed ICP profiles, and a handle to physically relate 
the ICP state to the DM's (see § 5.1). In the cases we 
have analyzed to now, we have identified two main clus- 
ter populations, HE and LE (see § 5.2). We have found 
the former to feature concentrations c « 3 — 5 associated 
with a slow outer decline of T(r) from a central plateau, 
a flat central and a steep outer brightness; one variant of 
this pattern is due to the ICP being rejuvenated by merg- 
ers, leading to the HE subclass. The other main popu- 
lation is constituted by the LE clusters. We have found 
these to feature: higher concentrations c > 6, associated 
with a central brightness spike and low but not vanishing 
central temperatures; and a steep outer decline of T{r) 
from the inner peak, with a considerable contribution of 
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turbulent support to equilibrium. Such a pattern is gen- 
erally sharpened toward low z (see § 6), and implies low 
outer SZ signals. Finally, our picture leads us to expect 
a main evolutionary sequence proceeding from HE to LE 
clusters. The above classification and time developments 
combine into our cluster Grand Design. 

In summary, in the articulated ensemble of galaxy 
clusters, the entropy-based framework provided by our 
Grand Design offers a thread toward understanding their 
basic astrophysics. Specifically, from the X-ray vantage 
point we interpret the correlations between ICP and DM 
parameters in terms of synchronization of the central 
and outer entropy demises, over timescales of several 
Gyrs. Within such a context, variance may be intro- 
duced by diverse large-scale environments adjacent to the 
outskirts, and possibly by multi-phase conditions at the 
center. Such a variance may blur the synchronized devel- 
opments, and originate instances intermediate between 
our two maijLclasses. We have identified one such en- 
semble, the HE clusters observed at relatively low z. At 
the other end z > 0.5, our Grand Design raises a specific 
issue concerning any clusters where central cooling is al- 
ready advanced (possibly requiring a multi-phase ICP), 
while entropy production is still high in briskly develop- 
ing outskirts (see also § 6). Observations of such objects 
at the current frontier of cluster astrophysics will consti- 
tute a challenging but rewarding aim. 
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APPENDIX 
CLUSTER BUILDUP 

In this Appendix we collect for the reader's convenience some basics of cluster formation, that are used throughout 
the main text to understand the entropy distribution in the ICP in connection with the halo development stages. Here 
our thrust will be to relate these stages to the shape of the initial density perturbation; on basic grounds, such times 

scale as tcoii {G p)^^/'^. 

In detail, the perturbation shape is conveniently parameterized as 5M /M cx M"*^, that may be considered a picccwise 
approximation to a realistically bell-shaped cold DM perturbation. Here 5M represents the mass excess within a shell 
at the initial comoving radius oc Af ^/'^ enclosing a mass M at background density. Such a shell will progressively 
detach from the Hubble flow, reach a maximimi 'turnaround' radius i?ta /{SM/M) oc and collapse back 

under local gravity to a virialization radius R « -Rta/2. 

The virialization occurs when SM/M attains the critical threshold 1.69 D~^{t) in terms of the linear growth factor 
D{t) depending on the cosmic time t. So the shape parameter e also governs the mass buildup after M{t) oc D^^'^{t) oc 
t'^/'^, where in the standard cosmology D{t) oc applies with d lowering from 2/3 to 1/2 as z decreases from above 
1 to below 0.5. The corresponding collapse time reads tcoii = M/M = et/d for the shell surrounding the mass M. 
Here e marks the cosmogonic effect of the perturbation tapering shape; on the other hand, d marks the effects of 
cosmology at large thinning out the background, and delaying collapse when d approaches 1/2. In many relations 
to follow what matters for the effective degree of halo development will be the combined index e/d = tcoii/t; values 
e/d < 1 apply to the fast collapse of the perturbation bulk, while during the slow outskirts development the accretion 
rate peters out corresponding to values e/d > 1. Note that the transition between the two regimes at Zt corresponds 
to e/d = 1, i.e., to the collapse time matching the Hubble expansion timescale as per definition; of course, this agrees 
with the transition epoch recognized in state-of-the-art A''— body simulations and semianalytic computations (see Zhao 
et al. 2003, Diemand et al. 2007, Fakhouri et al. 2010, Genel et al. 2010, Wang et al. 2011). At Zt, the halo 
concentration takes on the value c w 3.5 (with minor mass dependence, see Prada et al. 2011), and grows afterwards 
as c{z) 3.5 {1 + Zt)/ {I + z). 

While the halo develops and the concentration increases, entropy is produced by shocks in the contained ICP both 
at the center and at the boundary. As to the latter shocks, the driver is constituted by the external gas inflowing 
under the pull of the outer gravitational potential well; the inflow varies during cluster development, to the effects 
discussed below. 

Decreasing potential drops 
The potential drop from the turnaround i?ta to the shock radius Rs reads 

A$ = - / dr — — . (Al) 

2 ^ M'^-it „l-3e 



As the integrand behaves like 5M/r oc M '^/r oc r , one finds 

where Ac/i = A$/u^ is for the drop normalized to the circular velocity scale = GM(< R)/R, in fact at radius Rs 
of the boundary shock. The potential drop as a function of e is illustrated in Fig. Al. 

Outgrowing shock positions 
The position Rg of the shock may be determined from the scaling laws 

M 

vloz—A(p MocpviRl; (A3) 
Rs 

here M oc pR^ is the overall mass within Rs , Vi is the infall velocity in the cluster frame, p the background density 
(we have assumed ni (x p and m oc M), and At/) the adimensional potential drop described above. 
Combining the scaling laws yields 

i?. cx ^ (A<^)V3 « (A^)i/3i(^A+2)/3 (A4) 

that, when normalized to the turnaround radius i?ta oc M'^+^Z^^ j^ay be written in the form 

^oc(i)'^'(A<^)V3i-^+V3. (A5) 

as d takes on values within the narrow range 2/3 — 1/2, the explicit time dependence is very weak and may be neglected. 
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Using the expression of A0 derived in the previous § Al, the foUowing equation for x = Rs/Rta, obtains 

l_a;3e-2 3,_2 ' ^^^) 

the normahzation factor M is set on requiring that at the transition e = d the potential drop takes on the value 
A(j) « 0.57 corresponding to a 1.1 after Eq. (4) of the main text. This yields a shock radius i?s < i? « -Rta/2 close 
to the virial boundary during the early stages of cluster buildup that involve high accretion rates, corresponding to 
e<l. 

The position of the shock radius and the corresponding values of the potential drop are illustrated in Fig. Al. During 

the early collapse when e <C 1 applies wc find the approximations x oc e^^^ — > and A(/) cx e^^^^ to hold. At the 
other end, during the late outskirts development when e ^ 1 applies we obtain x — >■ 1 and Acj) rx e~^; thus the shock 
positions outgrows the virial boundary to approach the turnaround in the late development stage. 

Decreasing infall speeds and shock strengths 
Prom the scaling laws Eqs. (A3) we also derive an expression for the infall velocity (in the cluster frame) 

oc (a<^)V3 ^ (A<^)i/3 i(dA-i)/3 ^ (A7) 

This should be compared with the scaling = (5 fcsTi/S jimp)^^"^ cx p^^'^ of the sound speed in the preshock gas; from 
Eqs. (A3) we obtain p oc (M/M)"^ {A^)~^ for the density at the cluster edge, to yield 

)lV2/3 , , , —2/3 

c. ^ 1^ (A<^)- ^ (1) (A</,)-V3^-/3. (A8) 
The ratio of the two quantities reads 

(A<^)2/3t(<^A+i)/3 , (A9) 

and is seen to scale as £-4/5^d/3e e <c 1, and as e^^ t^^"^ for e ^ 1. In other words, strong shocks with vi ^ Cg 
take place during the early collapse of the cluster body, whilst during the late development of the outskirts the shocks 
weaken and vi <C Cg applies. 
From Eq. (A4) we compute the shock speed to be 

^ i?. oc (1)''' [I' + D (A<^)V3tW-i)/3 . (AlO) 
Taking the ratio of the two quantities (A7) and (AlO) yields the expression 

1)1 3 3 a 

This takes on values around 1/3 during the early collapse when e ^ 1; it grows during the late outskirts development 
when e ^ 1, since vi vanishes while Vg decreases toward its limiting value given by the sound speed Cg- Note the 
ubiquitous appearance in the DM dynamics of the key quantity A<f>, which will also appear directly in the ICP 
equilibrium condition. 

Weakening shocks and increasing seepage 

In the main text we discuss how the boundary shocks arc to weaken as the cluster outskirts develop; meanwhile, an 
increasing fraction of kinetic energy seeps through them. Next we explain why. 
The jump conditions for entropy, tempcratm-c and density across a shock front write (Landau & Lifshitz 1959) 

k2 ^ T2/Ti ^.^^ 
fci (n2/ni)2/3 

(A12) 

T2 _ 5 7 3 c2 n2 _ wi _ 4 

7\ ~ 16 ^ ^ 8 " 16 ' 1 + 3c2/52 ' 

The suffix 1 and 2 indicate pre- and postshock values, while quantities with a tilde refer to the shock rest frame (where 
the shock velocity Vs is zero by construction) ; in addition, = wi /c^ > 1 is the Mach number of the accretion shock. 
The behavior of these quantities as a function of A4 is illustrated in Fig. A2. 

In the cluster frame, the shock velocity Vs differs from zero, and the upstream and downstream bulk velocities are 
given by 5i,2 = ^^1,2 + Vs- Using Eqs. (AID) we work out the ratio V2/V1 to be 

vi 4 4 1 + Vs/vi 
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The above results are summarized as follows. During the early collapse of the cluster body with e <C 1, strong shocks 
with vi/cs » 1 and Vg — ui/3 hold; these imply high postshock temperatures fcsTb — M^^p oc e~^/^ and low bulk 
postshock velocities V2 oc e"^/^ ~ 0. On the other hand, during the late development of the cluster outskirts with e 1, 
weak shocks with Vi/cg <Si 1 and Vs — Cg occur, to yield low T2 ~ Ti and U2 — oc e~^. 

Thus as the cluster biiilclup progresses from bulk collapse to outskirts development, at the boundary the thermal 
postshock energy kBT2 monotonically decreases, but the bulk energy w| seeping through the shock to drive turbulence 
first increases up to a maximum; eventually, however, it decreases when the accretion becomes transonic. 



Specific clusters 

In the way of numerical values, the halo of a typical HE cluster observed at 2; « 0.1 collapsed at Zt ~ 0.5, developing 
a concentration c w 4. On the other hand, the halo of a typical LE cluster observed at « 0.1 collapsed at zt « 1; 
during the evolution e increases from initial values close to 0.3 to values around 0.6 at 2; w 0.3, and to 1 on moving 
to z « 0.1; meanwhile, its concentration increases from initial values around 3.5 to values around 5 at 2; w 0.3, and 
around 7 at z ~ 0.1. 

In parallel, for a relaxed, long-lived LE cluster the prevailing Mach numbers are to decline from « 10 to 3 and 
then toward 1 at low z, while the ratio increases from 10% to 30% with a correspondingly high outer turbulence 

levels, and then decreases again toward a few percent. 

Finally, one borderline instance is provided by the LE cluster A1689 at z 0.18; its high concentration c « 10 
implies the collapse redshift Zt « 2.5, particularly high for its mass M « 1.3 x 10^^ Mq. The other borderline instance 
is provided by the HE cluster A2218 again at 2: w 0.18; its concentration c w 5 implies the collapse redshift zt ~ 0.7. 

THE ICP SUPERMODEL 

Here we recap the basics of our ICP Supermodel introduced in § 3 (see Cavaliere et al. 2009, Fusco-Femiano et 
al. 2009, Lapi et al. 2010). The robust snapshots it provides guide our classification in cluster classes, and establish 
relationships betwcxai them to constitute the cluster Grand Design. 

In equilibrium conditions, the DM gravitational pull is withstood by the gradient of the thermal ICP pressure p, to 
yield 



GM{< r) 



1 dp{r) 
rUp n{r) dr 



5k^/^{r) _d 
2 /imp dr 



ksTir) 



(Bl) 



To begin with, in the second equality we have considered only thermal pressure p = n ksT/ /z, expressed in terms of the 
specific entropy k = ksT lv?l^\ in the main text the entropy distribution k{r) is related to definite physical processes: 
it is conserved by adiabatic compressions, produced at shock fronts, eroded by radiative cooling (nonthermal, turbulent 
support is dealt with in § B3). 

Eq. (Bl) shows how, given the potential well, the ICP disposition is set by the entropy distribution fc(r). In fact, it 
constitutes a first order linear differential equation for r(r), which solves to the profile 



T(r) = k^l^{f) 



1 



)fc-3/^(f') 



(B2) 



while the self-consistent density profile follows n(r) = [T{r)/k{r)]^^^. Here variables with a bar are normalized to their 
boundary values at r = R, while 6_r = iirup v^/ksTi expresses in the solution the boundary condition T2 provided by 
the shock jumps. 

The squared circular velocity v^{f) = M{< f)/f depends on the DM mass distribution. For the latter, we use our 
a-profiles derived from the Jeans equation (see Lapi & Cavaliere 2009a, 2009b, 2011); the corresponding density profile 
reads 



p(r) = r 



1 + w [cr)^ 



(B3) 



whore T « 0.76, 77 « 0.58, ^ w 4.56, and w = —(2 — r)/(2 — r — 77^) « 0.88 are constants with values suitable for rich 
clusters. Note that the standard NFW profile corresponds instead tor = l,7? = 0,^ = 2, and w = 1. It is seen how 
the density profile is modulated by c (the standard concentration parameter of the DM halo defined in § 1 of the main 
text), and particularly so in the outskirts. IDL and FORTRAN algorithms to implement the above equations can be 
foiuid at the URL http : //people . sissa. it/^lapi/Supermodel/. 

We stress that the standard models for the ICP distribution constitute useful approximations to the Supermodel 
results (e.g., Cavaliere & Fusco-Femiano 1976, Sarazin 1988); specifically, the isothermal /3— model applies in central 
regions with high central level kc, while the polytropic model k{r) oc n^~^^^{r) with index F w 1.2 applies in outer 
regions where n{r) drops quickly and the temperature T{r) undergoes a mild decline. 

In a nutshell, low entropy levels throughout the clusters allow the ICP thermal velocity dispersion T(r) to passively 
mirror the profile of DM velocity dispersion's cr^(r), as to share its radial run with a midle peak and a decline on 
both sides (as shown by Cavaliere et al. 2009, and confirmed by Hansen et al. 2010). This is because the DM and the 
passive ICP (in the absence of entropy additions) settle to a comparable equilibrium within the common gravitational 
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potential well. Conversely, high entropy levels at the center cause the ICP to resist gravitational compressions, and 
T(r) to maximally depart from a'^(r), so as to feature a monotonic increase inwards to a central plateau. 

Entropy distributions implemented 

Here we describe the radial entropy distributions k{r) implemented in the Supermodel, following the physical moti- 
vations given in § 2 of the main text. As also discussed there, which of the distributions is convenient to try first may 

be decided a priori from a quick look to the inner temperatures observed in X rays. 
For HE clusters (whose inner temperature profile is flat) we adopt the distribution 

fc(f) = fc, + (1 - fc,) f« ; (B4) 

this renders an outward rise with uniform slope a from a central level kc (see Fig. 1, red solid line). 
For HE clusters (whose inner temperature profile shows sharp wiggles) we adopt the distribution 

{kc for f < ff, 

(B5) 
fcc + (1 - kc) [(f - fy)/(l - ffT for f>ff . 

This represents a floor with level kc extending out to a radius ry, followed by an outward rise with uniform slope a 
(see Fig. 1, orange dashed line). 

For LE clusters (whose inner temperature profiles show an outward rise up to a maximum and then a decline) we 
implement the distribution 

{h {f/nY for f < n, 

, , (B6) 

where fcb = f^^° e" i'^-^b)/rb applies to ensure continuity at r = r^- This renders an outwards rise with uniform 
slope a out to a radius rh, and then a progressive bending follows with a linear decline of the slope, with the gradient 
a' = {a — an) fh/{—fh + 1) down to the boundary value an (see Fig. 1, blue dotted line). 

Observables, and parameter counting 

In connection with § 4, we detail how from T(r) and n(r) we proceed to compute the profiles of the X-ray and SZ 
observables. Specifically, the emission-weighted temperature is given by 

{T{w))=T2 ^° [l^ , (B7) 

Jo n2 A[T] 

where f = \/l — is expressed in terms of the projected radius w; here A[T] is the cooling function, with a typical 
dependence A oc T^/^ for hot clusters with average fegT > 2 keV. 
The brightness distribution is given by 

S{w)=nlR [ di n"^ A[f]F[EH,EL,T] ; (B8) 
Jo 

where the factor F[Eh, EL,T]c:i ^-E^/kBTir) _ ^-EH/kBT(r) ^^^^^^ ^^^^ account specific instrumental bands Eh - El- 
Fitting these expressions to the observations leads to pin down the following parameters and scales. From the profile 
normalizations we determine the ICP scales n2 and T^, and the DM scale R (if not independently given by observations 

of galaxy dynamics, 'red sequence' termination, gravitational lensing). 

From the profile shapes one can determine not only the ICP parameters describing the entropy run, but also the 
DM concentration parameter c, when not independently provided by gravitational lensing. For HE clusters 2 ICP 
parameters are needed, i.e., the central level kc and the slope a; for LE clusters 3 ICP parameters are needed, i.e., 

the outer value of the slope cr, the average derivative a' of the slope, and the radius r^; finally, for the HE clusters 3 
parameters are needed, i.e., the level kc and the extension of the central floor, and the outer slope a. 

A preliminary guideline as to which entropy shape is conveniently tried first is provided by a quick look at the gross 
temperature run at the center and in the outskirts, as discussed in § 4. A posteriori, the values of the reduced of 
the fits provide a final check. 

An independent observable is provided by the Comptonization parameter that marks the strength of the SZ eff'ect; 
it can be expressed as 

y{w) = n2T2R I di nf . (B9) 
Jo 

In the near future interferometric instrumentations like ALMA (sec http : / / science . nrao . edu/ alma/index . shtml) 
will provide measurements of comparable sensitive and resolution to the present X-ray instrumentations of XMM- 
Newton and Chandra. 
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Turbulent support 

In connection with § 3.2, we explain how the Supermodel can be readily extended to cover the ICP equilibrium 

when nonthermal, turbulent siipport contributes adding to thermal pressure (see Cavaliere et al. 2011). The relevant 
quantity is provided by the ratio S{r) = Pnth/pth of turbulent to thermal pressure or, equivalently, by the ratio 6/(l + S) 
of turbulent to total pressure ptoi = Pih (1 + <^)- 

We expect onset of turbulence in the outskirts of relaxed LE clusters (cf. Eq. A13), where weakening accretion 
shocks leave over an appreciable bulk energy to drive turbulent motions with maximal amplitude Sr w {v2/vi)^ up to 
30 — 40% at the virial radius. In fact, these motions start up with a coherence length L ~ R/2 set by the pressure scale 
height or by shock segmentation enforced by the adjoining filamentary structure, and then fragment into a dispersive 
cascade over the 'inertial range' to sizes i where dissipation begins. In the ICP context the dissipation scale writes 
^ ~ {c2/v)^^^ App (L/App)^/^ in terms of the ion coUisional mean free path App « 50 kpc and of the ratio v/c2 of the 
turbulent rms speed to the sound's. For subsonic turbulence with v/c2 < 1/3, the relevant scale £ exceeds somewhat 
App 100 kpc. This behavior may be rendered in terms of the simple functional shape 

<5(r)=^He-(«-'-)V^\ (BIO) 

which decays on the scale £ inward of a round plateau, a smoothed out representation of the inertial range. This 
provides the gradient of the turbulent pressure. 

In fact, the equation of hydrostatic equilibrium in the presence of turbulent support will contain the total pressure 
in the general form pth{r) [1 + S{r)], while the thermal component is still expressed as Pth{r) oc k{r) n^/^{r). On noting 
that ptot = Pth (1 + (5) = n ksT/ ji with T = T (1 + 5), it is convenient to introduce the extended entropy 

k = k{l + 5). (Bll) 

This quantity renders the conversion of kinetic energy into random energy at two levels, the microscopic one given by 
the standard k, and the dispersion into turbulent 'eddies' given by k 5. It is easily checked that in these terms the 
solution has the same form of Eq. (B2). 

It turns out that the profiles of emission-weighted temperature are little affected by turbulence (see Cavaliere et 
al. 2011); on the other hand, the thermal SZ effect is lowered to y = y/{l + 5). Finally, including turbulence in the 
equation of hydrostatic equilibrium brings the total mass reconstructed from X rays into agreement with the findings 
from simulations and with that measured via gravitational lensing observations (see Nagai et al. 2007, Lau et al. 2009, 
Meneghetti et al. 2010, Cavaliere et al. 2011). 
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Figure 1. The schematics illustrates our fiducial patterns for the ICP entropy distribution k{r). In the basic pattern (central level plus 
ramp; red solid line), entropy is raised at the boundary from intergalactic values fci ~ 10^ keV cm^ to high outer levels ^2 ~ 5 X 10^ keV 
cm-^ by strong boundary shocks. As the outskirts develop, the shocks weaken and the outer level lowers to k2 < 10^ keV cm^; meanwhile, 
the central entropy is eroded by radiative cooling down to low levels kc ~ 10^ keV cm^ (blue dotted line). On the other hand, blastwaves 
driven by deep mergers may enhance the central levels up to fcc ~ 3 X 10"^ keV cm^, spread out in the form of an extended floor (orange 
dashed line). 



16 



A. Cavaliere et al. 



A1 656 (Coma) z=0.02 



A2597 z=0.08 





A2256 z=0.06 



A1795 z=0.06 




0.1 

I- hi lUft} 



S 1 



steep LE \. 













i 











0,0 0,2 



0,5 0,8 1 




100 
r (arcsec) 



Figure 2. Supermodel fits to the briglitness and temperature profiles for the clusters A1656 (top left panel), A2597 (top right), A2256 
(bottom left), and A1795 (bottom right). Details are provided in Cavaliere et al. (2009), Fusco-Femiano et al. (2009), and Lapi et al. 
(2010). In the temperature panel of A1795, the Supermodel fit with a bending entropy profile is reported as a solid line, while the fit with 
a powcrlaw entropy profile is reported with a dashed line. 
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Entropy distributions 




r/R 

Figure 3. Entropy profiles normalized at -R500 as determined with the Supermodel for the 12 clusters listed in Table 1; red lines refer to 
HE and blue lines to LE clusters. In the radial range r > 0.2 R where cooling is negligible, these are overplotted to the outcomes of the 
nonradiative hydro-simulations by Burns et al. (2010); the black solid line represents the average over their sample of 24 relaxed clusters, 
with related variance illustrated by the shaded area. The four clusters presented in Fig. 2 are labeled. 
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Correlations 
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Figure 4. Top panel: central entropy level kc vs. entropy slope a in the cluster bulk. Dots illustrate our results from the Supermodel 
analysis of the 12 clusters listed in Table 1 (red dots refer to HEs and blue dots to LEs); squares are from the sample of 29 relaxed clusters 
by Pratt et al. (2010). Middle panel: central entropy level kc vs. outer entropy slope an; symbols are as above. Bottom panel: the DM 
concentration c vs. the outer ICP entropy slope an; symbols are as above. In all panels the Spearman's rank correlation coefficients p for 
the average data values are reported. 
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Figure 5. Outer temperature profiles T(r) we predict with the Supermodel from the existing inner/middle data concerning the clusters 
A1656 (left panel) and A2256 (right panel); data for the former are from Snowden et al. (2008) and Wik ct al. (2009), and for the latter 
from Snowden et al. (2008). The solid line represents our Supermodel fit in the region covered by the data, while the dashed line illustrates 
our prediction in the outskirts (for A1656 the fit has been performed basing solely on the inner data by Snowden et al. 2008). 
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Figure 6. The schematics illustrates the relationships among the cluster classes after the Grand Design. 
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Figure Al. Shock position /?s/_Rta relative to the turnaround radius and potential drop A(/> normalized to as a function of the 

parameter e; the inset illustrates the mass growth M(t) QC f^/" for three values of e = 1/6, 2/3, and 3, that span the range icoll/* = e/d 
going from 1/4 to 6. 
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Figure A2. Plot of the shock jumps A;2/fci, fl/j^j and T2/T1 as a function of the squared Mach number in the shock restframe; note 
that the divide between strong and weak shocks (associated with onset of turbulence) is around Ai^ = 3. During a cluster's evolution, the 
outskirts condition progresses from right to left. 
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Table 1 

A cluster library 



Cluster Class 


kc 


Tf/R 




n/R 


c 






[keV cm2] 


[xlO-2] 




[xlO-2] 







A1795t 

PKS0745t 

A2204t 

A1413t 

A2597t 

A2199t 

A1689t t 



LE 
LE 
LE 
LE 
LE 
LE 
LE 



151^ 

+6 
6 
+49 
49 



13 

105 



-3.43] 
-1.78^ 
0.16^ 
0.36"' 



3.36 
2.68 



^-2.68 
K+1.41 
-1.41 
K+0.31 
-0.31 
71 -H0.05 
"•'^-0.06 

n qc+o.oi 
<J-yo_o.oi 

n-H0.06 
-0.06 



28; 

23 
31 



-H2 
2 

-1-3 
3 
+ 7 
7 



27tl 



0.80"f 



8.5' 
7.6 
5.5 
8.3 
7.2 
6.7 
12.4 



1.9 
1.9 
1.7 
1.7 
1.1 
1.1 
1.7 
1.7 
+5.0 
5.2 
1.0 
1.0 
+5.3 
5.3 



0.3 (2.6) 
1.4 (4.4) 

1.1 (2.1) 

1.2 (1.9) 
0.3 
3.1 
1.7 



A2218* 


HE 


350111° 


A399* 


HE 


330ll«« 


A1656* 


HE 


,-20+160 
^-'"-160 


A644t 


iffi 


124+120 


A2256t * 


HE 


2481?^^ 



3 
12 



+0.2 
0.2 

+4 
4 



0.8; 
i.o; 

1.3 



+0.1 
0.1 
+0.1 
0.2 
+0.5 
0.2 



1 

-■■■^-0.1 



1.5 



+0.4 
0.3 



1+0.2 
-0.2 
1+1.5 
-1.5 
-,+0.8 
-0.8 
3+0.2 
-0.2 

2.7+1.7 



5.1^ 
3.1^ 
3.0j 
3.9l 



0.15 

1.3 

0.7 
0.5 (3.1) 
0.9 (2.5) 



Note. — Supermodel fits to the X-ray observables performed and/or refined in this work (marked with 
a *), and in the references Pusco-Femiano et al. (2009; marked with a X), and Lapi et al. (2010; marked 
with a t) . Dashes in the columns oirb/R {rf/R) indicate large (small) , irrelevant values. The last column 
provides the values of the reduced for the temperature fits, including r^- or r^, when necessary (in 
parentheses the values obtained on using the simple power law entropy run of Eq. 2). Note that for A2218, 
A399, A1656 the values rf = 1211, ^^0 1 ' ^-0 2' I'^spectively, come from the centrally flat brightness 
profile (see discussion by Fusco-Femiano et al. 2009). 



